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ABSTRACT 

The approximate nature of the adaptive-optics time constant measurements with MASS is examined. The 
cahbration coefficient C derived from numerical simulations of polychromatic scintillation shows dependence on 
. the height of the turbulence layer, wind speed, and seeing. The previously recommended value of C=1.27 is a 
good match to typical conditions, while C can vary from 0.6 to 1.6 in other circumstances. For two nights, MASS 
^ . was compared with the time constant measured with adaptive optics, and the expected agreement was found. 

We show that the single-layer approximation used in some AO systems to derive the AO time constant can give 
wrong results. A better approach is to estimate it from the speed of focus variation (the FADE method). The 
analysis of the speed of scintillation developed recently by V. Kornilov will lead to more accurate measurements 
of the AO time constant with MASS. 
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; 1. DEFINITIONS AND CONTEXT 

^ ', The adaptive optics atmospheric time constant tq depends on the wind speed and turbulence profiles. It is 
O ' defined as 
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where Cj^(h) is the vertical profile of the refractive index structure constant, V{h) is the vertical proffie of the 
modulus of the wind speed, h is the altitude above site. In the following we always assume that tq refers to the 
wavelength Aq = 0.5 /.im and observations at zenith. As the correction for the zenith distance depends on the 
unknown wind direction, it is simply ignored.'^ 

' The MASS instrument implements an approximate method® of estimating tq from the differential- exposure 

I scintillation index, DESI. The o'^esi computed for the smallest 2-cm MASS aperture as a differential index 
between 1 ms and 3 ms exposures. A formula 

TMASS - 0.175ms (a^Bs,)-"-' (2) 

has been suggested on the basis of limited data on real turbulence profiles. This formula is implemented in the 
, , standard MASS data processing 

V-1 ' found later by means of simulations that tmass calculated from ([5]) needs a corrective coefficient around 

I C = 1.27 The true time constant could be derived from the MASS data by applying this correction and adding 
the contribution of the ground layer (GL) which is not sensed by MASS (but measured with MASS-DIMM). 
Therefore, the AO time constant is estimated by MASS-DIMM as 

r-^/^ = {1.27 TMASS)-'^' + (0.057)-5/3Ao2 V^i' (C^d^GL- (3) 

The intrinsic accuracy of such estimate was evaluated to be ±20% or better. 

Travouillon et al^ have derived a somewhat different correction factor C = 1.73 by calculating tq from the 
turbulence profile measured by MASS and using NCEP wind velocities. Even larger factors of 2.45 and 2.11 
were determined in Ref. T'by the same method. This prompted a re-investigation of this matter by doing new 
simulations and comparing with alternative techniques. 
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2. NUMERICAL SIMULATIONS 



We simulate one phase screen at a given height h and calculate the intensity at the ground by means of the 
program simatmpoly . pro. The intensity pattern is a weighted sum for several wavelengths. Here we approximate 
the spectral response of MASS by 4 wavelengths of [400,450,500,550] nm with weights [0.31, 0.885, 0.60, 0.27]. 
This should mimic the response of TMT MASS-DIMMs, as studied by Kornilov.l^l He found for these instruments 
effective wavelength 474 nm and the bandwidth 99 nm (response curve without . crv). For our 4- wavelength 
approximation, the effective wavelength and bandwidth are 470 nm and 116 nm, respectively. The simulation 
does not include the light-source spectrum, assuming it flat. 



Figure 1. Fragment of the simulated scintillation pattern produced on the ground by a turbulent layer at 10 km with 0.5" 
seeing. The scintillation index is 0.27, the size of the fragment is 0.5 m. 

The simulated intensity distributions (Fig. [1} are saved in a binary file. In the previous simulator they were 
used in a Monte-Carlo approach where the 2-cm MASS aperture was "dragged" through the intensity screen, 
while 1-ms and 3-ms exposure time was emulated by suitable blurring of the aperture in one direction and by 
3x binning. Here we take a more direct approach and calculate DESI with the spatial filter 



P{fx) = smc{fxb) - sinc(3/j;6), 



(4) 



where fx is the component of spatial frequency, b — Vt^xp is the blur in the x-direction caused by the wind speed 
V during exposure time t^xp — 1 ms and sinc(a;) — sin(7ra::) / (ttx). The circular aperture of diameter d implies the 
filter 

TTfd 

The energy spectrum of the intensity is multiplied by the combined filter {AP)'^ and summed over all fre- 
quencies to get the DESI index. By omitting all filters, we obtain the raw scintillation index, by using only filter 
A - the scintillation index in the 2-cm aperture. Results for a test case were compared with previous simulations 
and found to be in agreement, validating the code. 

Considering that the calculation of intensity distribution is the most time-consuming task, we simulate the 
intensity screen for a given seeing and propagation distance and then calculate DESI and t^ass for a set of 12 
wind speeds, from lOm/s to 65m/s, by changing only the filter P. The calculation is repeated for 3 distances 
to the layer, 5, 10, and 15 km, and for 5 values of seeing e, from 0.3" to 1.5". Therefore we cover a wide range 
of conditions, with a varying degree of saturated scintillation. The largest scintillation index is 0.94 (layer at 



15kni, seeing 1.5"). In each case, the true time constant at Aq is known, tq = 0.31ro/y = 0.31(0.101/e)/y. We 
determine the corrective coefficient C — tq/tmass for each case. 





Figure 3. Left: Correction factor C as a function of the wind speed for a single turbulent layer at distance 15 km. Right: 
all values of tmass plotted against tq, the dotted line corresponds to C = 1.27. 



This study clearly shows that estimates of tq from DESI are approximate and can be off by as much as 
two times. No single corrective coefficient C can be derived, it varies from 0.6 to 1.6. For typical conditions 
{V = 30m/s, turbulence at 200mb), C — 1.27 appears to be a good choice. 



3. COMPARISON OF MASS WITH FADE 

Any working AO system provides real-time data on turbulence from which the AO time constant can be derived. 
For example, the covariance of atmospheric defocus falls to 1/2 of its maximum for a time lag to. 5 = Q.30D/V. 
It follows that To — l.Q5to,5{ro/D). A method of estimating tq from half-time correlation of Zernike aberrations 
has been proposed by Fusco et al.'^' and implemented in the NAOS at Paranal. It is valid for a single turbulent 
layer, but the authors argue that for several layers some "average" wind speed will result from to.5- This is not 
quite true, as we will see in a moment. 

TheorjPl'^ shows that the temporal structure function (SF) of defocus D4{t) produced by a combination of 



several atmospheric layers has the form 



Di{t) « 1.94i?5/3 ro;f 'X4(21/,t/D, e) + TV, (6) 

i 

where D is the telescope diameter, e is the central obscuration ratio, ro,i and Vi are the Fried parameters and 
wind speeds of the layers, N is the term caused by the measurement noise. The function K(x,y) is defined 
in Rcf. [5] It has a quadratic initial part at t ^ to.5i where small focus changes are proportional to time. In 
this short- exposure regime, the speed of defocus variation is related to the integral of V'^{h)Cf^{h)dh - the 2-nd 
moment of the wind speed V2- It is rather close to the 5/3rd moment V5/3 used in the definition of tq. The 
method to estimate tq from the speed of focus variation is called FADE (FAst DEfocus) It is resumed by a 
simple formula 

Di{t) « 0.036(i?/ro)-i/3(i/To)'. (7) 
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Figure 4. Left: example of the measured defocus SF (line) and its model by two layers (asterisks). The vertical axis 
m rad^ at 500 nm wavelength. Right: SF of defocus for two layers with 1" and 0.4" seeing moving at 10 and 60 m/s, 
respectively, the dotted line shows half-saturation and corresponds to to.s- 

Figure [3] shows the SF of defocus measured with the SOAR Adaptive Module (SAM) on October 2009. In 
this case D — 4.1m, e = 0.24. Taking the wind speed of 50 m/s, we estimate to.5 ^ 25 ms, or 6 loop cycles 
of SAM. Therefore, the temporal sampling of SAM (4.2 ms) is adequate for applying the FADE technique. By 
fitting the initial part of the SF with a two-layer model ([6]), we derive the defocus speed and tq. 

A simple test case presented in Fig. |4l right, shows why the half-time method can give wrong results. 
The to,5 = 0.105 s is essentially determined here by the strongest and slowest layer, leading to the estimate 
To = 2.39 ms. The actual time constant tq — 1.15 ms is mostly produced by the weak and fast layer which makes 
the dominant contribution to the speed of defocus and to the shape of the initial (quadratic) SF. The error of 
the half-time method to estimate To in this case is 2.1 times. 

So far, the SAM data were collected and processed for two nights, August 31 and October 2, 2009 (hereafter 
nights 1 and 2). AO loop data were recorded several times during ~ 12s. Atmospheric variations of defocus 
were derived from the signals of the corrector (DM), accounting for the frequency response of the closed AO 
loop. Voltages are transformed into Zernike coefficients in radians at Aq. The variance of the low-order Zernike 
coefficients is used to estimate the Fried parameter tq (seeing). 

Derivation of tq from the loop data involves few subtleties. First, we suppress all frequencies above 45 Hz 
because the data show focus vibration at 46 Hz (this low-pass filtering has little effect on the results, however). 
Second, the estimates depend on the maximum time lag used in the model-fitting. This length was varied from 
just 4 first points to 0.1s and 0.4 s. With increasing lag, the estimates tq^sam also become larger, see Table [TJ 
The fitted models show that on night 1 there was a fast-moving layer. 



Table 1. Average parameters from the FADE model fits 



Parameter 




Night 1 






Night 2 






4pt 


0.1s 


0.3s 


4pt 


0.1s 


0.3s 


to^mass/to^sam 


0.73 


0.83 


0.91 


1.42 


1.25 


1.02 


V, m/s 


30.5 


36.6 


27.6 


24.6 


27.8 


19.0 


Vlayer2, m/s 




63.9 


47.1 




64.9 


32.8 



The MASS-DIMM monitor at Cerro Pachon provided data for comparison with FADE. We used the ground 
wind speed and applied the formula (j3|) with C = 1.27 to get to_mass- The seeing on the night 1 was fast 
and bad, while the conditions on the night 2 (Fig. [SJ right) were closer to typical, but only 4 hours were clear. 
Different structure of turbulence on these two nights leads to different systematics between MASS and FADE 
(Fig. [51 left). We cannot tell whether C — 1.27 has to be increased or decreased, it is "about right". 



MASS-DIMM, Oct. 2/3, 2009 
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Figure 5. Left: Comparison of simultaneous ro from FADE (fitting 0.1s) and MASS; the dashed line shows equality. 
Right: Summary plot of atmospheric conditions on October 2, 2009 from the Cerro Pachon MASS-DIMM site monitor. 

4. CONCLUSIONS AND FURTHER WORK 

This study reveals, once again, the approximate nature of tq measurements from DESI implemented in the MASS 
software. Depending on the turbulence and wind profiles, the bias of MASS tq with the recommended correction 
factor C = 1.27 can be either positive or negative. Attempts to find the best calibration coefficient are hampered 
by the above-mentioned dependence on conditions. Our simulations show that the correction factor C = 1.73 
found by Travouillon et al.^^ is too high. However, the simulations do not account for the spectrum of the star; 
if it moves the effective spectral response of MASS blue-wards, the correction factor will be larger. 

Despite obvious shortcomings of the DESI method, it has a strong appeal, being a simple and no-cost addition 
to the MASS instrument. Considerable data on tq have been accumulated to date for many sites worldwide. 
We can do a better job on tq with MASS. Recent analysis by Kornilo\El uses the short-exposure approximation, 
where the dependence on wind is quadratic and separates from the spatial dependence. This permits to find 
optimum linear combinations of signals from MASS apertures and their pairwise covariances to cancel the height 
dependence, similarly to what is done for the free-atmosphere seeing and isoplanatic angle. By comparing this 
new method with DESI for one site, Kornilov finds C w 1.7. 

The second important improvement consists in using a longer effective exposure time when the turbulence is 
slow. This is possible with the new MASS software which records relevant statistical information. Unfortunately, 



this information is lost in the archival MASS data that cannot be cured post-factum; on good nights with slow 
turbulence when the DESI signal is weak, the tq derived from it is noisy and biasedP 

Finally, it becomes clear that for all turbulence-dependent optical parameters (scintillation, Zernike aber- 
rations, and even differential image motion) the speed of their temporal variation is directly related to the 
wind-speed 2nd moment ¥2, and hence can be used to measure tq. It is a practical matter to choose one or the 
other optical tracer for which the speed of variation can be measured accurately. The advantage of defocus in 
this respect is that, for a given aperture, it produces the strongest signal while remaining isotropic and immune 
to the instrument shake. But a DIMM with two or, preferably, 4 apertures could also be a promising solution 
for measuring tq, if a suitable theoretical analysis of its response is done. 
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